Anomalous Hall effects and electron polarizability 
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A theory of the anomalous and spin Hall effects, based on the space distribution of the current 
densities, is presented. Spin-orbit coupling gives rise to a space separation of the mass centers, as 
well as a current density separation of the quasiparticle states having opposite group velocities. It 
is shown that this microscopic property is essential for existence of both Hall effects. 
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It has been known for more than a century that a fer- 
romagnetic material exhibits, in addition to the standard 
Hall effect when placed in a magnetic field, an extraordi- 
nary Hall effect which does not vanish at zero magnetic 
field. The theory of this so-called anomalous Hall effect 
(AHE) has a long and confusing history, with different 
approaches giving in some cases conflicting results. While 
more recent calculations have somewhat unified the dif- 
ferent approaches and clarified the situation, it is still an 
active topic of research (see Ref. [f[ for a recent review). 
Closely related to the AHE is the spin Hall effect, rep- 
resented by spin accumulation on the edges of a current 
carrying sample @, 0] . 

It is generally accepted that the anomalous and spin 
Hall effects are induced by spin-orbit coupling. It was 
first suggested by Karplus and Luttinger [J] in 1954 to 
explain anomalous Hall effect observed on ferromagnetic 
crystals. Their analysis leads to the scattering indepen- 
dent off-diagonal components of the conductivity, which 
are assigned to the so-called "intrinsic" effect. Later, the- 
ories of this effect based on several specific models have 
been developed @, @ . As has been recently shown it is 
accompanied by strong orbital Hall effect [7, 8] . The con- 
ductivity is also affected by the scattering, which in the 
presence of spin-orbit coupling gives rise to so called side- 
jump [9] and skew scattering I3~12|. These also lead 
to anomalous Hall effect, called "extrinsic". The best 
quantitative agreement with experimental observations 
has been obtained with semi-classical transport theory 
[HI, leading to the Berry phase correction to the group 
velocity. For Fe crystals 14j it gives an anomalous con- 
ductivity ~ 750 fl^ 1 cm -1 while a value approaching 
1000 ft^ 1 cm -1 has been observed. 

It is the goal of this letter to shed light on AHE using a 
novel point of view. Our approach is based on the anal- 
ysis of the space distribution of local current densities, 
and it is simple and rather intuitive. We will show that 
the anomalous Hall conductivity is related to the spatial 
separation of the mass centers of states with opposite 
velocities. This confirms the interpretation of AHE in 
ferromagnetic systems as a consequence of the periodic 
field of electric dipoles (electric polarizability) induced 



by the applied current 0, EH EH , despite the fact that 
the original arguments were not convincing Q. We will 
also show that in non-magnetic systems, the spin-orbit 
coupling leads to a periodic variation of the spin polar- 
izability of the current densities in the transport regime. 
This effect can be viewed as an internal spin Hall effect. 
For the sake of simplicity, we limit our consideration to 
crystalline structures invariant under space inversion. 

Current density distribution is closely related to the 
orbital magnetization of solids. In zero magnetic field it 
has its origin in the orbital magnetic moment of atomic 
states. Ignoring spin effects, atomic wave functions in 
spherical polar-coordinate system can be written as 



yW(r,e,4>) = f a (r,6) 
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where m is the so called magnetic quantum number, with 
m — 0, ±1, • ■ ■ , and \m\ < I. It determines the magnetic 
moment along the z direction 
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where e denotes absolute value of the electron charge 
and v stands for velocity operator. The last expression 
represents a classical analogy with j a being the current 
flowing on a circular loop of the radius R a . Because of 
the energy degeneracy in m the total orbital magnetic 
moment vanishes. However, spin-orbit coupling together 
with exchange interaction remove this degeneracy giving 
rise to non-zero magnetic moment. 

Within a mean field approach the electron properties 
are controlled by a single electron Hamiltonian H con- 
taining two additive terms H so and H z representing spin- 
orbit coupling and an effective Zeeman-like spin splitting 
due to the exchange interaction, respectively: 



H = 



P 
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+ V(f) + H so + H z , 
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with mo being free electron mass, V(f) denotes the crys- 
talline potential, p is momentum operator and 
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VV(f)xp , H s = -fi s B eS -a, (4) 
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where A c denotes an effective Compton iength, and ele- 
ments of the vector a are Pauli matrices. Strength of the 
Zeeman-like splitting is controlled by the product of the 
Bohr magneton and the parameter B c g representing 
an effective magnetic field. The corresponding velocity 
operator reads 

" = — + tI^x VVYf) ■ (5) 
m Ah 

Eigenfunctions are spinors with two components, and 
since spin-orbit coupling does not destroy translation 
symmetry they are of the Bloch form. Energy spectrum 
E a (k) is a function of the wave vector k, with a being 
a band index, now including also, in addition to atomic 
orbital numbers, a spin number. Eigenfunctions are of 
the following form 

\a,k) = ^ a% {f) = ^Z=u a {k,¥), (6) 

and velocity expectation values are 

v a {k) = ^V i E a {k). (7) 

Spinors u a (k, r) are periodic functions of the lattice 
translation vectors. Assumed invariance under space in- 
version results in following fc-space symmetry: E a (k) = 
E a {—k) and v a (k) — —v a (—k). 

In order to analyze the role of the space distribution 
of the current densities, it is illustrative to present first 
the results for a simple model of a linear chain of atomic 
orbitals. It is assumed that this chain is forming a one- 
dimensional lattice along x direction with a period a. 
Model parameter will be chosen to satisfy conditions for 
which the tight-binding approach is applicable. Energy 
bands originated in overlap of atomic states |a) will be 
denoted by the corresponding magnetic quantum num- 
ber m. To model a ferromagnetic state, we assume that 
the effective field B c g is parallel with z direction and it is 
strong enough to ensure spin orientation along ±1 direc- 
tion with s z — ±1/2 being good quantum numbers. We 
have obtained numerical results by diagonalizing the sin- 
gle particle Hamiltonian for a two-dimensional separable 
chain potential 

V(x,y) = -V cos(2irx/a) + m nly 2 /2 . (8) 

Note that adding a z-dependent potential would not af- 
fect the current distribution of the considered model sys- 
tem. The parameters have been chosen to be in the tight- 
binding regime, i.e. to fully separate the studied band 
from the other energy bands (we used 2mQa 2 Vo/h 2 = 
75.0, m Q nla 2 = 1.4 x 4tt 2 V and irX 2 c /a 2 = 0.015). 

Typical current density distribution within the unit 
cell for the Bloch states \a,k x ) and \a,—k x ) are shown 
in Figs, la and lb, respectively. One observes circulating 
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FIG. 1: (color online) Current distributions for an energy 
band given by a chain of atomic states with spin s z = 1/2 
and m = — 1 for k x = 1.5/a (a,c) and k x = — 1.5/a (b,d). In a 
and b, the arrows indicate the direction of the current, and a 
lighter background indicates a larger current. In c and d the 
averaged current densities j"{k x ,y) are shown. 

currents forming vortices, which have the same orienta- 
tion for both cases. This orientation coincides with the 
orientation of the circulating current of the atomic or- 
bitals. In addition to this circulating current, there is a 
direct current flow, with opposite sign for the two cases. 
Corresponding total current and its direction is just de- 
termined by the velocity expectation value, —ev x (k x ). 
Note that these current flows are spatially separated in 
the two cases (they are on opposite sides of the circulat- 
ing current). The current densities averaged over x and 
z coordinates, 



jx( k x,y) = - 




clearly demonstrate the above mentioned spatial separa- 
tion of the currents having opposite velocity directions, 
as shown on Figs, lc and Id. 

This space separation of currents flowing in opposite 
directions is closely related to the mass-center separa- 
tion AY a (k x ) of states \ot, k x ) and \a, —k x ) with k x being 
consider as positive, k x > 0. For the considered model 
potential, Eq. ||5J), the y-component of the force operator 
reads 

1 A 2 
F y = —[p yi H] = -m ttly + 2s z -^m Cllp x . (10) 

In a stationary state the force expectation values has to 
vanish. Using the relation between p x and v x given by 
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Eq. ([5]) we get 



AY a (k x ) = 2s, , 
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In the limiting case of vanishing spin-orbit coupling 
AY a (k x ) — > and sum of the current densities j x (k x , y) + 
j x (~k x , y) approaches zero as well. 

Non-zero total current appears if there is different oc- 
cupation of states with opposite velocities which can be 
characterized by the chemical potential difference Afi. 
It can be related to the electric field along y direction, 
£ y = Afi/(e\AY a (k F )\), with AY a (k F ) being the mass- 
center separation of quasiparticles having opposite veloc- 
ities at the Fermi energy E F — /i. Within linear response 
approach the resulting current at zero temperature reads 
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Afi = -—AY a (k F )£ y . 

(12) 

Because of the non-zero separation AY a (k F ) and non- 
equal occupation of states with opposite velocities the 
applied current is giving rise an electric dipole moment, 
i.e. a charge polarization is induced. 

For the later use, let us express current J x (/i) in terms 
of the following quantity 



Pgjk) 



(a, k\r\a, k) x v a (k) , (13) 



where f2 ws defines a unit cell volume. Evaluation of the 
following expression 
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immediately gives the above result, Eq. (|12p . since ve- 
locity expectation value along y direction vanishes. The 
above defined quantity, Eq. (Q~3|), is the part of the or- 
bital magnetic moment within each of the unit cells which 
gives rise an electric dipole moment in the current carry- 
ing regime. For this reason it will be called as the orbital 
polarization moment. 

Generalization of the above treatment to a three- 
dimensional system is straightforward. Velocity expec- 
tation values have non-zero component also along y di- 
rection and they contribute to the orbital polarization 
moment defined by Eq. (|13p . The resulting contribution 
of the band a to the Hall conductivity component a xy 
can thus be written as follows 



ec 



5{E a (k)-n) P a (k) d 3 k, (15) 
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where integration is limited to the Brillouin zone and tt ws 
now denotes volume of the Wiegner-Seitz cell. Inserting 



for P a (k) and v a (k) their explicit forms, Eq. (TT3l) and 
Eq. ([7]), respectively, and using equality 
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already derived by Karplus and Luttinger [4], the inte- 
gration per parts gives the well known expression for the 
Hall conductivity of Bloch electrons 

2 p 

°-*»M = -7 f afcE f°( E «$)) \ti«(k)} d 3 k. (17) 
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Here fo(E) stands for zero-temperature Fermi-Dirac dis- 
tribution function and the Berry phase curvature fi a (fc) 
defined by the periodic part of Bloch functions, u a (k, r), 
reads 



Im ( Vc« a (fc,r) | x | Vr:u a (fc,r) 



(18) 



Our description of the anomalous Hall effect based on 
charge polarization effect is thus equivalent to the ap- 
proach based on the Berry phase correction fl3| | . 

For completeness, we must note that for the three- 
dimensional case, computing the quantity Ay, represent- 
ing the mass-center separation of states having opposite 
velocities along x direction, is not as easy as it was for 
the single atomic chain. It requires to express eigenfunc- 
tions in a mixed representation, to preserve the Bloch 
form along the x direction, while using Wannier repre- 
sentation along perpendicular directions. The Wannier 
representation gives functions which are bounded along 
the y direction, allowing to compute Ay. Although the 
explicit form of Ay is not simple, the main features are 
qualitatively the same as those presented for the single 
chain. 

Conductivity can directly be measured on samples hav- 
ing Corbino disc geometry. This arrangement can be 
modelled by considering a strip with periodic bound- 
ary conditions along x direction and current contacts at- 
tached to the strip edges allowing to apply current along 
y direction. It induces gradient of the electro-chemical 
potential along y direction and consequently it gives rise 
to the local charge polarization. As a result, Hall current 
is induced. Intra band scattering ensuring finite conduc- 
tivity across the strip is naturally of the side-jump char- 
acter. This type of scattering does not affect Hall con- 
ductivity which is given as the sum of the additive band 
contributions defined by Eq. (|15[) or Eq. (fP7j) . However, 
inter band scattering which is generally of the skew char- 
acter can affect Hall conductivity significantly. These 
features of the anomalous conductivity directly follows 
from the analysis based on quantum theory of the lin- 
ear response, Kubo formula. Although the derivation is 
straightforward, it is quite lengthy and will be thus pre- 
sented in a separate publication. 
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The relation we have presented here between the 
charge polarization and the anomalous Hall effect is sim- 
ilar to the discussion of Hall conductivity of Bloch elec- 
trons in rational quantizing magneticgv fields in terms of 
charge polarization. 17] However, in that case the phys- 
ical picture is strongly affected by chiral magnetic edge 
states. 
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FIG. 2: Transport current densities j^(k x ,y)Afj, for the en- 
ergy band of the two-fold degeneracy (B e s = 0) given by a 
chain of atomic states \a) with m = —1, s z = 1/2 (full line) 
and m — 1, s z — —1/2 (dashed line) for Fermi energy given 
by the wave number k x = 1.5/a. 



transport current density will show a periodic oscillation. 
This property can be interpreted as an internal spin Hall 
effect. At the sample edges of semiconductor systems, 
the oscillations of the spin polarizability will be modified 
by the confining potential defining sample edges. It can 
be expected that this modification is responsible for the 
already observed spin Hall effect. [U, Q 

To conclude, we have shown that the mass-center 
separation as well as the current density separation of 
states having opposite velocities is the essential feature 
of the systems with spin-orbit coupling. In the transport 
regime it gives rise to the charge polarization inducing 
anomalous Hall effect in ferromagnetic crystals. In non- 
magnetic systems it leads to a periodic spatial variation 
of the spin polarizability of the transport current density, 
predicted internal spin Hall effect, which is expected to 
be the origin of the spin accumulation at the edges of 
current currying samples. 
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Of particular interest are non-magnetic systems in 
which spin-orbit coupling is not negligible but effective 
Zeeman-like spin splitting vanishes, B c s — > 0. In this 
case the states of the single atomic chain \a,k x ) with 
orbital number m and spin s z are of the same energy 
as states \a, k x ) with opposite sign of the orbital number 
and spin, — m and — s z . Their orbital magnetization have 
opposite sign, the sum of their current densities vanishes, 

fx ( k x , y)+j € xi k x ,y)+j%(-k x , y) +j£(-k x , y) = 0, and to- 
tal magnetization vanishes as well. The mass-center sep- 
aration has also opposite signs, AY a (k x ) = —AY s (k x ), 
and in accord with Eq. (fT2|) the resulting anomalous Hall 
effect vanishes. However, the spin-orbit coupling has still 
an important effect in the transport regime. The current 
applied along the x direction gives rise for each band to 
non-equal occupation of states with opposite velocities 
represented by a local chemical potential difference A/x. 
The two considered bands have the same total current, 
but they have different space distribution because of the 
different mass-center positions determined by their spin 
orientation (Eq. (fTTj) ). As a result, the spin polarization 
of the transport current density averaged over x and z co- 
ordinates will be a function of the y coordinate. This is il- 
lustrated on Fig. [21 where the averaged transport current 
densities are shown for the same model parameters as in 
Fig. [1] Qualitatively the same features are expected for 
three-dimensional crystals: the spin polarization of the 
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